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Although the emission of radiation from dark matter annihilation is expected to be maximized 
at the Galactic Center, geometric factors and the presence of point-like and diffuse backgrounds 
make the choice of the angular window size to optimize the chance of a signal detection a non-trivial 
problem. We find that the best strategy is to focus on an annulus around the Galactic Center of 
~ 1° to > 30°, where the optimal size depends on the angular distribution of the signal and the 
backgrounds. Although our conclusions are general, we illustrate this point in the particular case 
of annihilation into two monochromatic photons in the phenomenologically most interesting range 
of energy 45 GeV < E < 80 GeV, which is of great interest for the GLAST satellite. We find for 
example that Dark Matter models with sufficiently strong line annihilation signals, like the Inert 
Doublet Model, may be detectable without or with reasonable boost factors. 
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I. INTRODUCTION 

Although the existence of Dark Matter seems 
to be required by a wealth of astrophysical and 
cosmological data (for reviews see [J 0), the 
identification of its physical nature remains still 
elusive. The most popular candidates are the 
so-called weakly interacting massive particles (or 
WIMPs), believed to be neutral and colorless ther- 
mal relics from the early universe with clcctroweak- 
scale masses (mx ~ 10—1000 GeV) and couplings. 
These generic assumptions seem indeed to roughly 
reproduce the amount of DM needed in the Uni- 
verse and in our Galaxy in particular, while arising 
at the same time naturally in many extensions of 
the Standard Model of Particle Physics. 

Among the different methods proposed to iden- 
tify WIMP dark matter (production in accelera- 
tors, direct detection via recoil of nuclei under- 
ground, indirect detection via their annihilation or 
decay products 0|), a particularly promising tech- 
nique is the detection of gamma rays from anni- 
hilation in the galactic halo. However, it is nowa- 
days clear that the DM signal should be likely sin- 
gled out of a larger astrophysical background, ei- 
ther point-like (the case of the Galactic Center) or 
diffuse. In assessing the capability of detecting a 
given WIMP candidate it is important to optimize 
the search strategy based on our current expec- 
tation for the energy range, but also the angular 
distribution of the signal. This information proves 
actually to be at least as useful in guiding DM 
searches as the energy spectrum and the choice 
of astrophysical targets. Here we illustrate this 
point with reference to the inner Galaxy region. 
Although the emission of radiation per unit solid 
angle from dark matter annihilation is expected to 
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be maximized at the Galactic Center, geometric 
factors and the presence of point-like and diffuse 
backgrounds make the choice of the window size 
to optimize the chance of a signal detection a non- 
trivial problem. To the best of our knowledge, this 
issue has been only raised within over-simplified as- 
sumptions for the signal and background [J or, by 
converse, treated in some technical studies by ex- 
perimental groups devoted to dark matter searches 
(see for example H), with the background profile 
determined e.g. via a numerical model. In no case 
a comparative study for different profiles and back- 
ground assumptions, which illustrates the physical 
origin of the optimal cuts to adopt, has been pre- 
sented. This is the main task we address in this 
article. Wc shall prove that — in all the commonly 
considered models for the DM distribution in our 
Galaxy and given the present knowledge of point- 
like and diffuse backgrounds — the optimal strategy 
for DM searches is never to focus on the inner 0.1°- 
0.5° around the Galactic Center. Unless astro- 
physical processes dramatically enhance the DM 
distribution in the innermost region of our Galaxy, 
producing effectively a bright point-source, an an- 
nulus of radius from ~ 1° to > 30° is always to be 
preferred. The optimal size depends on the angu- 
lar distribution of the signal and the backgrounds, 
and the improvement in the sensitivity reaches a 
factor of a few, thus affecting the prospects for the 
detectability of specific models. To give an explicit 
example, we shall refer mostly to the detection of 
gamma-ray lines from annihilation. Gamma-ray 
lines of energy E = mx would be an unambiguous 
signature of WIMPs as DM particles populating 
our Galactic Halo, but the two-photon channel is 
loop(or radiativcly)-supprcssed with respect to the 
main WIMP annihilation channels: fcrmion pairs, 
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weak gauge bosons or higgses. Even models where 
the line emission is enhanced, as the Inert Doublet 
Model (IDM) @, are very challenging to detect 
via line emission. Yet, with an optimized choice of 
the angular size, the IDM model may be detectable 
without or with reasonable boost factors (i.e., en- 
hancement of the flux with respect to the smooth 
halo signal due to the "dumpiness" of DM) . Simi- 
lar considerations apply for the more conventional 
continuum gamma-ray flux from DM annihilation, 
usually considered a more promising channel. 

This paper is structured as follows. In Sec. [U 
we discuss the ingredients contributing to the DM 
signal. In Sec. IIIII we describe our model for the 
point-like and diffuse backgrounds, while Sec. IIVI 
is devoted to present our results. In Sec. |V] we 
conclude. 



II. THE SIGNAL 

It is customary to write the number flux of 
gamma rays (per energy and solid angle) produced 
in dark matter (DM) annihilations in our Galactic 
halo as 
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where the factor in square brackets depends only 
on the particle physics model: (av) is the WIMP 
annihilation cross section multiplied by the rela- 
tive velocity of the two WIMPs (averaged over the 
WIMP velocity distribution), mx is the mass of 
the WIMP, dN^/dE is the photon spectrum per 
annihilation. If the dark matter is not its own 
antiparticle (as assumed here), Eq. (fl]) should be 
multiplied further by a factor 1/2 (if X and X are 
equally abundant). 

The factor outside the brackets only depends on 
the DM astrophysics, in particular on the distri- 
bution of DM in the halo and on the position of 
the Solar System within the halo. The quantities 
df and pt are constants with dimensions of dis- 
tance and energy density respectively, whose only 
purpose is of making J — which encodes the angu- 
lar dependence of the DM signal — a dimensionless 
function. Following Q , we fix these parameters at 
df = 8.5 kpc and pi = 0.3GeV/cm 3 (indicative the 
galactocentric distance of the sun and of the local 
density of dark matter). The J- function is thus 
defined in terms of the DM density p as 
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In the previous equations, ip is the angle between 
the direction in the sky and the Galactic Center 
so that, in terms of galactic latitude b and lon- 
gitude Z, cos-0 = cos b cos I; r© « 8.0 kpc 0] is 
the solar distance from the Galactic Center; I is 
the distance from the Sun along the line-of-sight 
(l.o.s.). Note that we implicitly assumed that the 
halo is spherically symmetric and that the angular 
dependence only comes from the off-center position 
of the Sun in the halo, which implies for example 
the azimuthal symmetry of the signal (see e.g. [lj| 
and refs. therein for a discussion of sub-leading 
effects determining the angular distribution of the 
signal). A general class of spherically symmetric, 
smooth halo distributions can be fitted as 
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where po is a normalization constant and r s is a 
characteristic radius below which the profile scales 
as r~ 7 . Two of the most well known profiles 
have been proposed by Navarro, Frenk and White 
(NFW) [ll| and Moore et al. [HJ. Another fit, 
proposed by Kravtsov et al. [H| mimics a flat core 
potential, which seems to be in better agreement 
with some observations. More recent simulations 
(see e.g. 0, EE [H, H3) suggest that halo density 
profiles are better represented by a function with a 
continuously-varying slope. At least above the re- 
solved radius (~ 1 kpc) and for Galactic-size halos, 
the newly proposed fitting formulae provide only 
marginal improvement with respect to the more 
traditional ones, although it seems that the inner 
slope of a DM halo is very likely shallower that 
the Moore one, and probably shallower than the 
NFW one, too. For illustrative purposes, we shall 
consider the Moore, NFW and Kravtsov profiles 
with the parameters reported in Table I, since the 
former and the latter bracket the different cases, 
with the NFW (as well as the paramterization used 
in [l7|) falling in between. Similarly to Ref. [H[, 
the normalizations are chosen so that the mass con- 
tained within the solar circle provides the appropri- 
ate DM contribution to the local rotational curves, 
such that pq = p(r®) — the DM density at the so- 
lar distance from the GC — yields to the canonical 
value p Q = 0.3GcV/cm 3 for the NFW profile. 



TABLE I: Halo parameters, see text and Ref. [TJ 



Model 


Q 
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PQ 


r s 










[GeVcm" 3 ] 


[kpc] 


Moore 


1.5 


3 


1.5 


0.27 


28 


NFW 


1.0 


3 


1.0 


0.30 


20 


Kravtsov 


2.0 


3 


0.4 


0.47 


10 



Note that in none of the above-mentioned mod- 
els astrophysical processes modifying the pure DM 
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profile are included. For example, "adiabatic com- 
pression" may lead to the formation of a DM spike 
around the central supermassive black hole; merg- 
ers as well as scattering on the dense stellar cusp 
around the central black hole may destroy density 
enhancements, etc. While potentially relevant in 
determining the dark matter annihilation bright- 
ness of the Galactic Center, they only affect re- 
gions too close to the Galactic Center to be re- 
solved angularly by present detectors. In general, 
they might lead to a measurable gamma ray flux 
from the innermost angular bin even in presence 
of relatively large astrophysical backgrounds, but 
in this case only the energy information is relevant 
for the separation, for which the present discussion 
is irrelevant. We not address these effects further; 
details can be found in the rcfs. cited in [l[ or 
ref. [3. 

The spectrum per annihilation depends on the 
nature of the WIMP, although the gamma-ray 
yield is often largely dominated by the result of 
the fragmentation and hadronization of the WIMP 
annihilation products. Neutralinos in the MSSM, 
for example, typically annihilate to final states 
consisting of heavy fermions and gauge or Higgs 



bosons 



With the exception of the r + r chan- 



nel, each of these annihilation modes result in a 
very similar spectrum of gamma rays. In most of 
the following considerations, however, it is the an- 
gular distribution of the DM radiation that deter- 
mine the optimal angular window. To avoid unnec- 
essary complications, we shall illustrate our point 
considering mostly the case of annihilation in two 
gammas, for which 
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where & 77 is the branching ratio for the 2-photon 
line. Although it would constitute a striking sig- 
nature of a DM annihilation process, in general 
the photon line is considered hard to detect, since 
it is loop-suppressed with respect to the main an- 
nihilation modes Q. This is not necessarily the 
case in some DM models, as for example the In- 
ert Doublet Model (IDM) 0,0. In general, the 
line feature may be enhanced whenever: (i) the 
DM almost decouples from fermions, while being 
mostly coupled to gauge bosons; (ii) the DM can- 
didate is in the mass range 45 GeV ~ mz/2 < 
Mx < mw — 80 GeV, the lower limit coming 
from the bound on the Z°-width measurement and 
the upper limit to avoid the annihilation mode 
XX —> W + W~ by kinematical constraints (other 
cases where a prominent line emission is possible 
at lower or higher energies have been treated re- 
spectively in 21 1 and [22| ). Interestingly, this is 
probably the best region where the GLAST satel- 
lite [23[ has a reasonable chance to detect a line 
emission and at the same time, potentially interest- 



ing for the next generation of low-threshold atmo- 
spheric cherenkov telescopes ( ACTs) , like CTA [24| 
or AGIS H. In the following, for GLAST we 
shall use an integrated area exposure over time of 
10 4 cm 2 yr, corresponding to about 5 years of ob- 
servations, and an energy resolution of ±7% at la. 
For ACTs, we shall assume an effective array area 
of A c s = 1 km 2 , energy resolution of 15% and 200 
hour observations, with a 2° x 2° field of view. 

Whenever needed, we shall compare our results 
with two benchmark models for IDM candidates 
among the ones considered in 0, see Table HI 
Also, to get a feeling of how things change if the 
continuum gamma-ray spectrum is considered in- 
stead, we shall also report the results obtained 
when a fiducial case of a 100 GeV MSSM neu- 
tralino mainly annihilating into W + W~ is consid- 
ered. 

TABLE II: Benchmark models in the IDM scenario, 
see 0] ■ 



Model 


mx 


(<rv) 






[GeV] 


[cm 3 /s] 




I 


70 


1.6 x 10~ 28 


0.36 


II 


50 


8.2 x 10" 29 


0.29 



III. THE BACKGROUNDS 

The Galactic Center is a complex region of the 
sky at all wavelengths, the gamma-ray window be- 
ing no exception. In this section, we discuss how, 
in our analysis, we treat the backgrounds for dark 
matter searches due to known and unknown astro- 
physical sources of gamma rays. Since we focus on 
an extended signal, diffuse (or unresolved) back- 
grounds are our main concern. However, to illus- 
trate how effective is a large field-of-view search, 
we shall compare it with the sensitivity achieved 
when focusing on the inner region around the 
Galactic Center, when at least one astrophysical 
point-like background is known to exist. This is 
the relatively bright, very high-energy gamma ray 
source observed by HESS, MAGIC, WHIPPLE, 
and CANGAROO-II (26|. This source is consis- 
tent with point-like emission and is located at I = 
359°56'41.1" ± 6.4" (stat.), b = -0°2'39.2" ± 5.9" 
(stat.) with a systematic pointing error of 28" (27l ], 
It appears to be coincident with the position of Sgr 
A*, the black hole constituting the dynamical cen- 
ter of the Milky Way. Following the measurements 
of HESS, we describe the spectrum of this source 
as a power-law given by: 

/ rp \ -2.25 

$ GC -1.0xl0- 8 f-f-l GeV" 1 cm- 2 s- 1 . 
\GeVJ 

(6) 
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At energies below ~ 200 GeV, the spectrum of this 
source has not yet been measured. GLAST, how- 
ever, will be capable of measuring the spectrum 
of this source at energies below the thresholds of 
HESS and other ACTs. 

In addition to the HESS source, a yet uniden- 
tified source has been detected by EGRET in 
the GeV range, approximately 0.2° away from 
the dynamical center of our galaxy [28|, |2{|. It 
has not been detected so-far at the energies ex- 
plored by HESS and other ACTs, which means 
that its power-law spectrum cuts-off between ~10 
and ~100 GeV. Although it should be included in 
more realistic analyses of the DM signal from the 
inner galactic regions [3(J, in the following we ne- 
glect it, thus overestimating the diagnostic power 
of a narrow field of view search and providing a 
conservative estimate of the improvement achieved 
by a large angle search. 



The overall diffuse gamma-ray radiation can be 
qualitatively divided into a galactic and an ex- 
tragalactic contribution. Since the latter is not 
simply the isotropic part of the flux, the sepa- 
ration of these two components can be done at 
present only assuming a specific model for the 
production of secondaries by cosmic rays in the 
galactic disk and halo. (However, a measurement 
of the cosmological Compton- Getting effect that 
should be achievable for GLAST would provide a 
model- independent way to separate the two contri- 
butions [3l|). A significant fraction of the quasi- 
isotropic component, especially in the GeV range, 
may be due to high-latitude galactic emission com- 
ing from processes in the magnetized halo of the 
Milky Way. Wc employ here a fit of the galac- 
tic diffuse flux proposed in Q and calibrated on 
EGRET data around the GeV [H, 



I gal (E) = N (l, b) x 10" 6 (7^) cm-W^GeV- 1 , 
where the arguments I, b are in degrees, -180° < I < 180° and -90° < b < 90°, a ~ 2.7, and 
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The EGRET collaboration derived the intensity of the extragalactic gamma-ray flux as [32 



I m {E) = fc x 10" 



E 



0.451GeV 



'- 2 s- 1 8r- 1 GeV _1 , 



(9) 



with 13 = 2.10 ± 0.03, k = (7.32 ± 0.34) and the fits is valid from E_~10 MeV to E ~ 100 GeV. 



In the following, we shall consider two models 
of diffuse backgrounds: a "conservative" model 
assuming a = 2.7 in Eq. ([7]) and ko = 7.32 in 
Eq. @ , and an "optimistic" model, where a harder 
spectrum a = 3.0 is assumed for the Galactic 
background and k^ = 3.66. Based on presently 
available models, one expects indeed that GLAST 
should resolve roughly half of the extragalactic 
background ; analogously, one expects that un- 
resolved sources contribute significantly to the dif- 
fuse Galactic spectrum as measured by EGRET, 
and once detected and removed by GLAST the 
Galactic background should be steeper (see e.g. 
discussion in [2l|). 

Finally, for ACTs the dominant background is 
the isotropic cosmic ray one. At energies well be- 
low ~ 100 GeV one expects the cosmic-ray electron 
background to play a dominant role. In the follow- 
ing, we account for these backgrounds following the 



formulae reported in 



IV. RESULTS 

First, let us repeat a simple argument already 
sketched in 0]. The DM signal coming from in- 
tegrating over a solid angle AQ is proportional to 
(J(Q,))ACl, where 



(10) 



For an annulus around the GC with inner open- 
ing semi-angle ip m and outer one ^m, the previous 
formula reduces to 
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FIG. 1: Plot of the function {J)ip M defined in Eq. JTTJ 
(proportional to the significance of a DM signal for 
an isotropic background in the limit S <C B) for the 
three models considered in Table I: the solid lines as- 
sume ipm — 0.1°, the dashed lines tpM — 0.3°; the left- 
most (red) lines are for the Moore profile, the rightmost 
(green) lines for the Kravtsov, in the middle (blue) lines 
for the NFW profile. 



The significance S/N of a Signal S in presence 
of a background B is defined as 
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The dark matter signal goes like S oc (J) AO, while 
for an isotropic background B oc AO. In the limit 
S -^i B (which is likely to hold in all realistic DM 
searches), we have 

-|oc<J>VAO~ (J)ip M . (13) 

In Fig. [TJ we plot this function vs. ipM for 
the three models considered in Table II, assum- 
ing ip m = 0.1° (solid lines) or xp m = 0.3° (dashed 
lines). Note the exact choice of the inner bound is 
irrelevant for cored profiles: the peak of this func- 
tion is above 20°. Even for profiles as cusped as 
the NFW (or more) the peak of this function is at 
a few degrees. Of course, the Moore profile gives 
a formally divergent signal when ip m — > 0, so the 
inner angle cutoff does matter. 

Figure [TJ already shows that the S/N is typi- 
cally maximized for values of ipM 3> V'm- However, 
the diffuse backgrounds are not {I, &}-indcpcndent, 
but really decline away from the Galactic Plane. 
When using the models for the galactic and ex- 
tragalactic background introduced in the previous 
section, we see how this effect is appreciable. In 
Fig. [U , in order to make the comparison among 



the different models easier, we plot the rescaled 
S/N vs. 6 max (in such a way that their maximum 
values are equal) for a region 0.4° < |6| < 6 max , 
0° < |/| < Z max = Lax for the fiducial IDM model 
I of Tab. II, and for the continuum spectrum 
of a 100 GeV MSSM ncutralino annihilating into 
W + W~. In the former case, we consider the 1-er 
smoothing of the gamma-line due to the GLAST 
energy resolution, while for the continuum we ac- 
tually show the S/N in the most significant energy 
bin (which happens to fall in the ~ 1 GeV range) . 
The inner cut has been chosen to 0.4° to avoid con- 
tamination from known point sources close to the 
GC. Note how the peak of S/N moves to larger an- 
gles in all cases compared with Fig. [TJ developing 
a maximum at finite angle also for the Moore pro- 
file. For the continuum spectrum the peak of the 
significance moves to larger angles, due to the fact 
that the isotropic component of the background is 
harder, and thus at the lower energies relevant for 
this signal the overall background is dominated by 
the Galactic background which falls rapidly with b. 
It is also worth nothing that even for heavy dark 
matter particles (TeV mass scale) the continuum 
photon spectrum peaks at or below about 30 GeV. 
Since the S/N for continuum spectra is dominated 
by its value at the peak, the angular window re- 
sults reported in Figs. 1 and 2 are really more 
general than the case they were discussed for, re- 
flecting the situation for the continuum spectrum 
from DM in all the interesting mass range. 

In Fig. [3] we show the curves with S/N = 3 for 
the branching ratio into two gammas, normalized 
in terms of the fiducial annihilation cross section 
value (av) = 3 x 10 -26 cm 3 /s in the phenomenolog- 
ically most interesting mass range below mw- The 
triangles represent the fiducial points of the IDM 
model reported in Table II, and for the window size 
that maximizes S/N (respectively & max ~ 1.5°, 15° 
and 30°, for Moore, NFW and Kravtsov profiles of 
Tab. I). The solid lines are for the conservative 
background estimate, the dashed lines for the op- 
timistic one. It is interesting to note that for the 
Moore profile, the fiducial models should be de- 
tectable at the 3(7 level without need for boost 
factors. For the NFW case, boost factors of order 
10 should suffice. Even for the unfortunate case 
of a Kravtsov profile, a boost factor ~ 20 should 
be sufficient. These values ap pea r to be viable in 
cold DM cosmologies, see e.g. [3^]. Thus, our con- 
clusions are more promising that what envisioned 
in 0], where boost factors > 100 were considered. 
Also, the top dot-dashed line shows the sensitivity 
achieved if one only includes the inner 0.1° region 
around the GC, where the background is domi- 
nated by the central source. These limits, com- 
parable with the results reported in [36j | . show an 
evident improvement. 

It is also worth commenting on the fact that, 
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FIG. 2: relative Signal/Noise vs. annulus size i> max (we 
equate their maximum value to ease the comparison) 
for a region 0.4° < |6| < 6 max , 0° < |Z| < Z max = 
bmax for the continuum spectrum of a 100 GeV MSSM 
neutralino annihilating into W + W~ (dot-dashed line, 
NFW profile) and for the fiducial IDM model I of Tab. 
II for the Moore, NFW and Kravtsov profile (solid lines 
respectively from the left to the right). 

according to our estimate, although ACTs are not 
optimized to search for line signals in the range be- 
low 100 GeV, the performance of the next gener- 
ation, km 2 instruments presently being considered 
may be comparable or better than GLAST. For ex- 
ample, in accordance with existing estimates (37j . 
for the fiducial exposure parameters we have used 
and the NFW profile we find a sensitivity a fac- 
tor ~T.5 better in the parameter space of Fig. 3, 
despite their worse energy resolution. It is also 
fortunate that the ground-based instruments, that 
have necessarily a much smaller field of view with 
respect to satellite observatories like GLAST, do 
not suffer too much from limited field of view, at 
least for relatively cusped profiles. The dominant 
backgrounds of ACTs are indeed isotropically dis- 
tributed, and for this case the optimal annulus size 
is < 2°-^3° (for cusped profiles), within the field of 
view range considered by present designs for future 
ACT. 



V. DISCUSSION AND CONCLUSIONS 

A large field of view instrument like GLAST 
promises to revolutionize the field of high energy 
gamma ray astrophysics. But it could also turn 
into a discovery instrument for DM, if it is made 
by WIMPs which annihilate into standard model 
particles in the halo of our galaxy. Since the DM 
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FIG. 3: Plot of the curves with S/N = 3 in the 10-100 
GeV range, with vertical lines in correspondence of 
mz/2 and mw- The dot-dashed line shows the sensi- 
tivity achieved if one only includes the inner 0.1° region 
around the GC for the NFW profile. The couples of 
solid/dashed lines (solid for conservative background, 
dashed for optimistic one) represent instead the sensi- 
tivity for the optimal annulus size for each model: from 
top to bottom, for Kravtsov, NFW and Moore profiles, 
respectively. 



signal grows like the square of its density, it is 
usually assumed that the Galactic Center, where 
the DM density is expected to be the highest, is 
the most promising region to focus on. However, 
two arguments conspire to modify this naive expec- 
tation: i) the presence of astrophysical point-like 
backgrounds close to the GC; ii) geometric con- 
siderations, in particular related to the different 
angular shape of the DM signal and astrophysical 
backgrounds. 

It is thus important to establish the level of 
the diffuse galactic and extragalactic diffuse back- 
grounds, a task which will be accomplished with 
exquisite precision by GLAST itself. Yet, by em- 
ploying simple models based on EGRET data, we 
showed that it is always to be preferred a window 
size > 1° around the GC, even for profiles as steep 
as the Moore one. The flatter the profile, the larger 
the optimal window size: for example we find an 
optimal galactic latitude cut fr max > 10° for NFW, 
frmax si 30° for Kravtsov profile. The discovery 
of additional point-like sources is not expected to 
change things dramatically, since especially at high 
energies the good angular resolution of GLAST 
will allow one to remove the quasi-totality of the 
resolved point-source contaminations with only a 
minor loss in the solid angle of the diffuse signal. 

When these geometric effects are taken into 
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account and the search strategy is adjusted con- 
sistently, the prospects for DM detection improve 
considerably. For example, we found that models 
with sufficiently strong line annihilation signals, 
like the Inert Doublet Model, may be detectable 
without or with still reasonable boost factors, of 
the order of 10 even for Kravtsov-likc profiles. 
Also, the perspectives to detect a DM signal 
become less sensitive to the unknown profile of the 
DM halo. In particular the S/N of quasi-cored 
profiles is only a factor ~ 2 worse than for the 
NFW case. We also showed with one example 
that our considerations remain true also for the 
widely considered case of continuum gamma-ray 
spectrum from DM annihilation. Given the 



generality of these considerations, we believe they 
should be taken into account when estimating the 
reach of DM searches via gamma-ray telescopes. 
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